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Abstract 

Consider the infinite rate mutually catalytic branching process (IMUB) constructed in [ID] and |11] , 
We show that for finite initial condition only one type survives in the long run if the interaction kernel is 
recurrent. On the other hand, under a slightly stronger condition than transience, we show that both types 
can coexist. 
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1 Introduction and Main Results 
1.1 Background and Motivation 

In [5] Dawson and Perkins considered the following system of coupled stochastic differential equations as a 
model for mutually catalytic branching: 

Y tit (k) = Y ifi (k) + [ J2 A ( k > l ) Y iAl)ds (f.f) 
Jo ies 



+ ( 7 yi, s (fc)y 2 , s (fc)) 1/2 dW itS (k) for t>0,*eS,i = l,2. 
Jo 

Here S is a countable set that is thought of as the site space and 7 > is a parameter. (In fact, Dawson and 
Perkins made the explicit choice S = Z d .) The matrix A is defined by 

A(k,l) = a{k,l)-l {k=l} , (1.2) 

where a is the transpose of a stochastic matrix a T on S such that sup fcG 5 X);es a (k, I) < 00. Note that A T is 
the g-matrix of the continuous time Markov chain on S with jump kernel a T . (In fact, Dawson and Perkins 
assumed that a is symmetric, but this is not substantial.) Finally, (W^(fc), k € S, i = 1, 2) is an independent 
family of one-dimensional Brownian motions. 



1 This work is partly funded by the German Israeli Foundation with grant number G-807-227. 6/2003 



1 



This is a spatial model for the evolution of two populations i = 1, 2. Yi :t (k) is the size of the population of type 
i at site k G S at time i. The individuals migrate on the site space S according to the discrete space heat flow 
induced by a. Furthermore, at each given site each type of the population undergoes a random dynamic that 
can be interpreted as continuous state Feller's branching with a branching rate proportional to the local size 
of the respective other type. Both types undergo (independently) the same branching dynamics and influence 
each other in a symmetric way - hence the name mutually catalytic branching process. 

Dawson and Perkins studied the longtime behaviour of this model with summable initial conditions and estab- 
lished a dichotomy of coexistence versus non-coexistence of types depending on transience and recurrence of the 
Markov chain associated with a. Via the self-duality of the mutually catalytic branching process, its total mass 
behaviour for summable initial conditions provides information on the local behaviour if the initial condition is 
infinite and sufficiently homogeneous. For x S [0,oo) 2 , let x denote the state in ([0,oo) 2 ) s with x^k) = Xi for 
all k E S, i = 1,2. Assume Yq — x. It is shown in [3 Theorem 1.4] that Y t converges in distribution to some 
random field Yoo as t — > oo. Furthermore (under some mild regularity assumptions on A T ), we have 

Px[>i,oo (0)^2,00 (0) > 0] > a T is transient. 

Hence in the recurrent case, the two types segregate locally and form clusters. The assumption that the initial 
point is constant can be weakened to an ergodic random initial condition (see [3]). 

The starting point for this work was the wish to get a quantitative description of the cluster growth in the 
recurrent case. We only briefly give the heuristics. Dawson and Perkins also constructed a version of their 
process in continuous space R instead of S as the solution of the stochastic partial differential equation 

dY- t(r) I ■ 

g ' = AF M (r) + ^<yY u (r)Y 2 , t (r)Wi(t,r) for r el, i= 1,2, (1.3) 

where W\ and W 2 are independent space time white noises and A is the Laplace operator. As A on R is 
recurrent, here also types segregate. Now due to Brownian scaling, if we denote by Y 1 the solution of (|1.3[) 
with that given value of 7, we obtain 

P^(rVT)) reR e .] = Pj(>7 T M) reR e •]. (1-4) 

Equation (|1.4p shows that clusters of Y^t grow like y/T and that a better understanding of the precise cluster 
formation can be obtained by letting 7 — > 00 for fixed time. We call the limiting process A of Y 1 as 7 — > 00 
the infinite rate mutually catalytic branching process (IMUB). 

In [9] the infinite rate mutually catalytic branching process X was constructed for 5 a singleton. It was shown 
that Y 1 converges to X as 7 — * 00 and X was characterised in terms of a martingale problem and in terms of 
its generator. Since the two types cannot coexist in the limit 7 — > 00, the proper state space for the one-colony 
IMUB is 

E := [0,oo) 2 \ (0,oo) 2 . 

In [TU] the IMUB process X was constructed for countable 5* via approximate solutions to a related Poisson 
noise stochastic partial differential equation. Here A is a strong Markov process that takes values in a suitable 
subspace of E s that fulfills some growth condition (Liggett-Spitzer space). More precisely, as shown in [TOl 
discussion after (1.32)], for k € S, i = 1, 2, the process 

Mi, t (fc) := X i>t (k) - X u0 (k) - f AX i>s (k)ds, f>0, (1.5) 

Jo 

is an LP martingale for every pe [1,2) (but not for p = 2) that could be represented as the stochastic integral 
with respect to a Poisson point process. 

Furthermore, the IMUB process A was characterised as the solution to a certain martingale problem and it was 
shown that Y 1 converges to A. In order to formulate the martingale problem, we need the notation 

xoy := -{x\ + x 2 ){yi + y%) + i{x\ - x 2 )(yi - V<z) for x,y 6 E 
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(with i — y— T) and 

for x, y G such that the sum is well-defined. Then the E s valued Markov process X with initial value 
Xq = x is characterised by the requirement that 

((AX s ,y))e« x '<y» da, t > 0, (1.6) 
be a martingale for all suitable y € S s . 

In [TT] a construction of X is performed via a Trotter type approximation scheme, see also [T3]. Loosely 
speaking, for given e > 0, consider a process X s that solves (jl.lj) with 7 = in each interval [tie, (n + l)e), 
n € No. At the times ne the state X^ e _ is replaced by the limit (as t — + 00) of a solution Y of (jl.ip with 
.4 = and 7 > with initial state Y := X^ e _. That is, the value X^ e _(k) at each colony k S S is replaced 
(independently) by a point in E chosen randomly according to the exit distribution of planar Brownian motion 
in [0, oo) 2 started at X^ e _(k). (See Section |2~21 for a more detailed description.) It was shown in [TT] that X e 
converges as e — ► to a process that solves the martingale problem (|1.6p . 

In this paper, we aim at understanding the longtime behaviour of X for summable initial states x (and thus, in 
(|1.6p we could take any bounded y € E s ). Let r be the amount of time that two independent Markov chains 
with ^-matrix A T spend together. We show that if r has infinite expectation (and A T fulfils some modest 
regularity condition), then the types cannot coexist in the long run. On the other hand, if A fulfills a condition 
that is somewhat stronger than E[r] < 00, then the types can coexist in the long run. As X is a process with 
an infinite variance random dynamics, the meta theorem that relates stability of the longtime behaviour to 
transience of the migration dynamics does not apply here. It remains open to check if there are cases where 
E[t] < 00 but coexistence of types is impossible. In particular, it would be interesting to know if this could 
happen for certain (transient) random walk kernels A. 



AXuv)) 



A*>y)) 



1.2 Results 

Before we present our result, we give a more detailed description of the result on the longtime behaviour in [5]. 
Let 

MT t := (Y ht , 1) := ]T Y i>t (k), t > 0, (1.7) 
kes 

be the total mass process of type i = 1,2 and assume that M^ ,M^ S (0, 00). Since and M% are 
(orthogonal) nonnegative martingales they converge almost surely and in fact also in L 1 . Denote by Mj x the 
limit variables. For the case where a is a random walk kernel, Theorem 1.2 of [5] takes the concise form 

p [ M i y oo M 2 Y oo > 0] >° a is transient. (1.8) 

In order to formulate the result for the more general case, we have to be a bit more careful. Let at be the 
continuous time kernel; that is 

°° in An 00 i.n n n 

a t = cMAt) :=E 1 T- = E e- tL ^-, (1-9) 
£ — ' n! L — ' n\ 

where a n and A n denote the matrix powers. Note that af is then the continuous time Markov kernel related 
to a T . Furthermore, for t > 0, define the Green kernels 

G t (k,l):= / a s (k,l)ds and G{k,l) := / a s (k,l)ds. (1.10) 
Jo Jo 
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Finally, let 

G* := sup G t (k,k). (1.11) 

fees 

We say that coexistence of types is possible if 'P[MY 00 M)f 00 > 0] > for all initial states Y with MY ,M 2Q G 
(0,oo). We say that coexistence of types is impossible if PiMYoo^Yoo > 0] = for all initial states Y with 
M£ ,M£ e[0,oo). 

Now Theorem 1.2 of [5] states the following: 

Theorem (Theorem 1.2 of [5j) Assume that a is symmetric, 
(i) If swp keS G(k, k) < oo, then coexistence of types is possible, 
(n) If 

inf liminf^ilMl >0 (1.12) 

fees t^oo Gl 

and if a is recurrent and irreducible, then coexistence of types is impossible. 

The theorem describes a dichotomy between a stable behaviour (coexistence of types) and an instable or clus- 
tering behaviour (segregation of types) depending on potential properties of the underlying migration dynamics. 
A similar dichotomy along the same line of transience and recurrence (in the case of migration of random walk 
type) was observed before for many interacting models with finite variance dynamics such as the voter model 
(see [TJ [7]), interacting diffusions on a compact interval ( [ISl EE H] ) , branching random walk ([8]), the gener- 
alized smoothing and potlatch process process ([6]), and the so-called linear systems ([12l Chapter IX]). When 
the local dynamics has moments only up to order 1 +f3 for some (3 £ (0, 1), then the the critical line between the 
two regimes may shift to the point where higher powers of the Green operator are finite (see [4]). In the model 
we study in this paper, all moments less than the second are finite and thus a dichotomy could be expected 
close to the transience/recurrence line but a little shifted to the transient side. 

In order to prepare for the formulation of our theorem and since we want to get rid of the symmetry assumption 
that Dawson and Perkins imposed on a, we have to introduce some more notation first. Recall that a T (k,l) = 
a(l, k) is the transpose of a and define af , Gf and so on similarly as for a. Furthermore, define the symmetrised 
kernels 

a t (k,l) = af /2 a t / 2 (k,l) = a t/2 (m,k)a t/2 (m,l) 

meS (1.13) 
- 1 1 
A=-(A + A T ), a=-(a + a T ). 

Note that (<it)t>o is a semigroup (and is generated by A) if AA T = A 1 A. In particular, if S is an Abelian 
group and a is a random walk kernel, then (a t ) t >o is the semigroup of the difference of two (rate 1/2) random 
walks and its one step transition matrix is a. 

Define the expected amount of time r two independent Markov chains with q-matrix A T spend together: 

G{k,l) := lim G t (k,l), where G t (k,l) := [ a s (k,l)ds. (1.14) 

As mentioned above, for many finite variance models, finiteness of G is equivalent to stability. We will show for 
the IMUB model here that G(0, 0) = oo (plus some mild regularity conditions) is sufficient for non-coexistence 
of types. In order to formulate the regularity conditions properly, we will also need 

G* t := supG t (fc, k). (1.15) 

fees 

For showing coexistence of types we need more refined quantities. Define 

a s (k, I) — (a s (a + a T )aJ) (fc, I) — a s (k,m)(a(m,n) + a(n,m))a s (l,n). (1-16) 

m.neS 
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Let us give a very rough heuristic, why this object appears here. If we start with unit mass of type 1 at 
ki and unit mass of type 2 at fe, then, as we will show, the expected mass of type i at time s at site m is 
a s (m, ki). Recall that the types exclude each other at any given site. If type i is absent at m at time s, then 
the infinitesimal impact of type i at site m is governed by the immigration at rate a(m, I) from the other sites 
I G S; that is, it is of order aa s (m, ki). Summing over all sites, we see that the expected total "activity" is of 
order a s (k\, k 2 ). Since the interaction of types has infinite variance, it is not a itself that is the crucial quantity, 
but rather we will see that we will need a logarithmic correction term. In order quantify the total amount of 
interaction in the "transient" case, we define 



G log (k,l):= a s (k,l)(l + \\og(a s (k,l))\)ds. (1.17) 



It is easy to check that 



Giog(M) > / a s {k,l)ds = G(k,l) - t {k=l} . 







Hence G\ og is infinite if G is infinite. 
Define the total mass process 

M itt :=(X itt ,l) (1.18) 

and assume that Mo G [0, oo) 2 . Recall the martingales Mj(fc) from (|1.5p and note that M.^ t = YlkeS ^,t(^) 
since YlkeS ■A(h,l) = for all I £ S. Hence Mi and M 2 are is a nonnegative martingales and therefore the 
almost sure limits 

Z % := lim Ma 

t — >oo 

exist, £ = 1,2. 

Recall that a T is a stochastic matrix and that A(k,l) = a{k,l) — l{k=i}- 

Theorem 1 Let X be the infinite rate mutually catalytic branching process with kernel a. Assume the initial 
state is 

Xi,o = !{/!} and X 2 ^ = t{i 2 } 

for some € S, h 7^/2- 

If l\ and I2 are such that G\ og (l\, I2) is sufficiently small, then there is coexistence of types in the longtime limit; 
that is, P[Z 1 > 0, Z 2 > 0] > 0. 

Theorem 2 Assume that a fulfills 

supa(/c,fc)<l (1.19) 

kes 

c:= inf liminf Gt( - fc,/) > 0. (1.20) 

Then coexistence of types is impossible; that is, if Mq € [0,oo) 2 , then Z\Z 2 = almost surely. 

Remark 3 (i) Note that in the case where a = a T , condition (|1.20p implies irreducibility of a. Furthermore, 
if a = a T is recurrent and irreducible, due to the strong Markov property, (|1.20[) is equivalent to (|1.12[) . 

(ii) Assume that S is an Abelian group and a is a random walk kernel. In this case (|1.20|) is equivalent to a 
being irreducible and recurrent. In particular, (11.201) holds if a is irreducible and recurrent. 

(iii) For symmetric simple random walk on the d-dimensional integer lattice Z d , d > 3, it is simple to show 
(using the CLT), that Gi og (ii, ^2) ~ Cd\\h — h\\ 2 ~ d as — Z2II — ► 00, for some constant Cd € (0, 00). 
Hence the assumptions of Theorem Q] are fulfilled for simple random walk with the two populations being 
sufficiently far apart. 



and 
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The proofs in [5] heavily rely on second moment methods. The main difficulty in the proofs here is the lack of 
second moments. (For this reason, presumably the statement of Theorem [1] fails under the weaker assumption 
that only G < oo.) The strategy of proof for Theorem [1] is therefore to introduce for K > an auxiliary process 
X K whose jumps in each coordinate are suppressed when they lead out of the square [0,-ftT] 2 . This is done in 
such a way that the coordinate processes (minus the drifts) become square integrable orthogonal martingales. 
For these martingales we use the conditions on Gi og to estimate the conditional quadratic variation process. 

The proof of Theorem [2] also uses the auxiliary process X K and its conditional quadratic variation process, but 
the arguments are more involved. 



1.3 Organization of the paper 

In Section[2]we prove Theorem[T] We start with deriving basic properties of Brownian motion in [0, K] 2 stopped 
upon hitting the boundary such as hitting distribution, moments and so on. Then we construct the auxiliary 
process X K and conclude Theorem Q] via second moment estimates. 

In Section [3] we prove Theorem [2j 



2 Coexistence of types, Proof of Theorem [T] 

In Section 12. II we make preliminary calculations on the variance of planar Brownian motion in [0,-ftf] 2 stopped 
upon hitting the boundary. In Section l2~2l we construct the auxiliary process X K . In Section l2~3l we show that 
the coordinates of X K are orthogonal martingales and we compute their conditional quadratic variation. In the 
final Section 12.41 we put the ends together to conclude the proof of Theorem [TJ 

2.1 Brownian Motion in a Square 

We start with the definition of some spaces. Recall that E = [0, oo) 2 \ (0, oo) 2 and define 

E K ■= En[0,K) 2 , E K ■= [0,K} 2 \(0,K) 2 , E K := E K \E K . (2.1) 



Denote by Q the harmonic measure of planar Brownian motion in [0, oo) 2 . That is, if B = {B\,B2) is a 
Brownian motion in R 2 started at x G [0, oo) 2 and 

T := in£{t > : B t ^(0,oo) 2 }, 



then we define 



(2.2) 

If x = (u, v) G (0, oo) 2 , then the harmonic measure Q x has a one-dimensional Lebesgue density on E that can 
be computed explicitly 



4 uvu 
n Au 2 v 2 + (u 2 + 
4 uvv 



v 2 — u 2 ) 



n Au 2 v 2 + (v 2 + u 2 - v 2 Y 



du, if v — 0, 



dv, if u = 0. 



(2.3) 



Furthermore, trivially we have Q x = S x if x G E. 
In [5] it is shown that for x G (0, oo) 2 and p G (0, 2), 



Vi Qx(dy) = Xi 



(2.4) 



G 



and 

E4r p/2 ]<oo. (2.5) 
Applying the Burkholder-Davis-Gundy inequality this implies that 

JylQ x {dy)<™. (2.6) 

However, for p = 2 we have 

r y?Qx(dy) =00. (2.7) 



Now consider planar Brownian motion B on [0, K] 2 and its harmonic measure Q K defined by 

Qf := P X [B TK g •], 

where 

r K :=mf{t>0: B t ^{0,K) 2 }. 
Due to the obvious scaling, we can restrict ourselves mostly to K = 1. For simplicity, we define r := t%. 

Lemma 2.1 For all x E [0, l] 2 and i = 1,2, rae Ziaue 

[y i Ql(dy) = x i . (2.8) 

Fitrt/iermore, 

F(x) := E a [r] = / ( W - Z,) 2 Qi(dy) (2.9) 



and 

'(i/i-si)(ife-a;2)Q£(<ii/) =0. (2.10) 



Proof. (-Bj )tAr ) t >o is a bounded martingale and f At is its quadratic variation. Hence, (|2.8[) and (|2.9p are 
simple consequences of the optional stopping theorem for martingales. Similarly, (|2.10p follows from the fact 
that the product (-Bi jtAr -B2,tAT)t>o i s a bounded martingale. □ 

Lemma 2.2 For x — (u,v) € [0, l] 2 , V(u,v) has the Fourier expansion 

oc oo 

E( U) „)[r] = = ^ ^a mi „ sin ((2m + 1)ttu) sin ((2n + 1)W) , (2.11) 

m— n— 

w/iere 



32 



a ™.« : = ^ 9 fTWXTfcrT^ f° r al1 m >» e N o- (2.12) 
7r 4 (2m + l) z + [2n + \) z 

Proof. It is well known that V is the unique solution of the Poisson equation with Dirichlct boundary condition 

^Ag= -1 in (0, l) 2 , 

2 y v ' ; ' (2.13) 

g = at 5(0, l) 2 . 

Denote by g(u, v) the right hand side in (|2.1ip . Clearly, g = at 9(0, l) 2 . Furthermore, for u, v S (0, 1), 

W-> - ill £ =^r^) <"<> 

\ m=0 / V n=0 / 

where the last equality follows from the fact that each factor is the Fourier series of the function on (0, 1) that 
is constant 1. Hence, we have g = V. □ 
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Corollary 2.3 For planar Brownian motion B in [0, AT] 2 , we have 

V [u , v) {t k }=K 2 V(u/K, v/K) 

and 

Cov (UtV) [B i>TK ,B jjTK ] = K 2 V(u/K, v/K) t {z=j} . 
Proof. This follows from Brownian scaling. □ 

Lemma 2.4 For all u, v £ [0, 1], we have 

V(u,v)>2u(l-u)v(l-v). (2.15) 

Proof. Let f(u, v) = 2u(l - u) v(l - v). Then 

-Af(u, v) = -2[u(l -u) + v(l - v)] > -1. 

Hence, by the maximum principle, / is a sub-solution for the Poisson problem (|2. 13|) which shows / < V. □ 

Proposition 2.5 Let K > and let B be Brownian motion in [0,K] 2 . Then we have 

Cov (u , v) [B i>TK , B j!TK ] <8««[l + log(lO + (Iog(l/t4)Alog(l/t;))] l {4=i} for all u,v£[0,K] (2.16) 

and 

Cov {u , v) [B i<TK , B jtTK ] > -uv for all u,v £ [0,K/2]. (2.17) 



Proof. It is enough to consider the case K = 1. Note that (I2.17|) is an immediate consequence of Lemma [2^41 
Hence we concentrate on showing (|2.1G|) . 

We have to show that (with V from (|2.11| ) 

V(u,v) < 8w[l + log(l/u)]. (2.18) 

By symmetry in u and v, this implies (|2.16p . 

As V(u, v) is bounded by the expected time one-dimensional Brownian motion started at v needs to hit {0, 1}, 
we have V(u,v) <v(l-v)<v for all u,v £ [0, 1]. Hence for u > 1/3, (|2~T5|) holds with C = 3/(l + log(3)) < 2. 

We may and will now assume u < 1/3. Let M £ N be such that 2A j +3 < u — 2M+1 ■ ^ e w ^ show that 

V(u,v) < 8w[l + log(M)] for v £ [0, 1], M £ N. (2.19) 

We estimate | sin((2m + l)7ru)| < 7r(2m + 1) u for m < M — 1 and | sin((2m + l)7ru)| < 1 for m > M, as well as 
| sin((2n + l)7rv)| < ir(2n + 1) v. Hence we obtain 

oo M-l 

V(u, v) = a m ,n sin ((2m + l)iru) sin ((2n + l)irv) 

n— m— 

oo oo 

+ a m .ri sin ((2m + l)7Ttt) sin ((2n + l)irv) 

n=0 m=M 

32uv%^^ 1 32v 

< _ 



Ex -» 1 32v x 1 x -> 1 

^ (2m + l) 2 + (2n+ l) 2 + T 5 " ^ 2m + 1 <W2m + l) 2 + (2n + l) 2 

m=0 n=0 v ' v ' m=M n=Q y ' v ' 

I M {u,v) + J M (u,v). 
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The two summands will be estimated separately. First note that 

M-l „oo 



, , 32m; ^ f°° 1 

1m{U, V) < TT~ > / — —5 -TT dt 

TT 2 ' J 2m +1 ) 2 + t 2 

m— 

16 ™ ^i/-^ 1 1 16 uv r 
= > < 1 + log M) . 

m— 

Similarly, we get (note that 5M > 2M + 3 > 1/u by the assumption on M) 

T . , 16 1; ^ 1 4„ 1 20w 1 20 

m=M v ' 

Summing up and noting that 16/7T + 20/ir 2 < 8, we obtain (|2.19[) . □ 



2.2 Construction of the Truncated Process 

The aim of this section is to construct a process X K that approaches X as K — > 00 and which has finite second 
moments. The idea is to suppress the large jumps of X so that the remaining jumps have second moments. It 
turns out that if we proceed a bit more subtly, then we can obtain even that the coordinate processes of X K are 
orthogonal square integrable martingales and that we can control the conditional quadratic variation process. 
The rough idea is as follows. The jumps of X can be interpreted as being driven by the positional changes of 
planar Brownian motion at its exit points from [0,oo) 2 . For the process X K , we stop this planar Brownian 
motion when it exits [0, if] 2 . 

We could proceed in two ways to construct X K : 

1. We could imitate the SPDE construction of X (see [10]) by replacing the intensity measure on E of the 
Poisson point process by a suitable intensity measure on [0, K] 2 \ (0, K) 2 . 

2. We could imitate the Trotter type construction of X (see |11| ) by replacing the harmonic measure Q on 
[0, oo) 2 by the harmonic measure Q K on [0, K] 2 . 

Here we follow the latter approach. In [TT] the following was done in order to construct X: For fixed e > 0, 
consider the stochastic process X 6 with values in ([0,oo) 2 ) s with the following dynamics: 

(i) Within each time interval [ne, (n + l)e), n € No, X 6 is the solution of 

dXl t (k) = (AXf t ){k) dt for t e [ne, (n + l)e) k £ S. 
Clearly, the explicit solution is 

Xt t (k) = (at- ns Xl nE )(k) for t £ [ne, (n + l)e). 

(ii) At time ne, X s has a discontinuity. Independently, each coordinate X £ ne _(k) = Oe-X"( n _j\ e (&) is replaced 
by a random element of E drawn according to the distribution Q x ■ (k) ■ 

In order that the solution in Step (i) is well-defined we have to impose some growth condition on the initial states 
(see [11] Theorem 1]). Since here we are interested in finite initial states, this growth condition is automatically 
fulfilled. 

In [11] it was shown that X s converges as e — > to X in the Skorohod space of paths [0, 00) — > ([0, oo) 2 ) s . 
Define 

t k ■= inf {t > : (X ltt + X 2>u 1) > K/2}. 
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Clearly, tk is a stopping time and by Doob's inequality, 

P[r K < ooj < 2 — . 

IS. 

Now assume that (x± +X2, 1) < if- We construct X ' £ with initial condition x just as X but with two differences: 

(1) In Step (ii) above we replace E by E K and Q by Q K . 

(2) If (X*£ e + X^ n E £ , 1) > K/2, then Step (ii) above is omitted. 

Note that Step (i) preserves the total mass, hence once the total mass exceeds K/2, the process X K,e is simply 
the discrete space A heat flow. Denote by 

r K ' s := inf {t > : (X*' e + X*' e , 1) > K/2} 

the time when this first happens. Note that due to the strong Markov property of planar Brownian motion, we 
have 

Q x = [ Q*{dy)Q y . 
Je k 

Hence X s and X K ' £ can be coupled to coincide almost surely until T K - e . Since X s converges, this implies that 
also (X K jf ) t >o converges in the Skorohod space as e — > 0. Since the A heat flow clearly exists, in fact X K,e 
converges as e — > to some process X" . Clearly, 

M%°(k) := X«f{k) - f (AX«f)(k) ds, i = 1,2, k e S, 
Jo 

are orthogonal square integrable martingales with conditional quadratic variation process 

n: ne<tAr K - E 

Upper bound for the conditional quadratic variation 

Note that for each n with [n — l)e < t K ' £ , we either have X 1 ^_i\ s {k) = (which implies X^ £ _(k) < 
eaX^ £ i _ 1 ^ £ (k)) or X% jt-i\ s {k) = (which implies X 2 ^_(k) < sa X 2 ^^(k)). Hence (denoting by I the 
identity matrix) 

n: ne<tAr K ' E 

where (by Proposition [53]) 

< h K (u) := 8Ku(l + log(K/u)). (2.20) 

Since bounded L 2 martingales converge to bounded L 2 martingales, and since s~ 1 (a £ — I)(k,l) a(k,l) for 
k 7^ I, we get that 

M£(fc) := X$(k) - / (AX?)(k)ds, i — 1,2, k <E S, (2.21) 
Jo 

are orthogonal square integrable martingales with conditional quadratic variation processes 

(Mf(k)) t < f (aX« s (k)h K (X* s (k)) + aX^ s {k)h K {X^ s {k))) ds, i = l,2, fceS. (2.22) 
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Lower bound for the conditional quadratic variation 

Define 

t k := inf {t > : {Xg t + X% t , 1) > K/2}. (2.23) 

By Doob's inequality, 

v[t k <oo] < 2 (j^g + ^2,o,i) j 

iv 

Furthermore, X coincides with X (in distribution) until time r . By Proposition [231 we have 

n: ne<tAr K > e 

Recall that a £ (kj) = e~ e + e~ £ ea(k 1 1) + • . Since ^ (^_ 1 v(fc)^2 < (^_ 1 ) £ .(fc) = f° r an n < t X ' £ , we have 

<Mf-(fc)) t > lee-' £ [^S-D^W^-DeW + ^-i)^) ^-^WL 

n: ne<fAT K ' £ 

where we also used aX^' E n (fc)Xf ' £ , .,(&) = n (fc)xf' £ , n (&) for fn - l)e < t K ' £ . 

Since (X^' £ ) e> o is a convergent sequence of bounded square integrable martingales, also the conditional 
quadratic variation processes converge and we infer for t > s > 0, 

{Mf(k)) t - (Mf (k)) s > - / (aXK r (k)XK r (k) + aX« r (k)X« r (k)) dr. (2.24) 

L Js/\T K V ' 

2.3 Truncated Process and Martingales 

In order not to interrupt the flow of the argument later, we start here with a lemma. 

Lemma 2.6 Let Y and Z be non- positively correlated non-negative random variables and assume that h : 
[0, oo ) — ► [0, oo) is concave and monotone increasing. Then E[Yh(Z)] < E\Y] /i(E[Z]). 

Proof. If E[Z] = 0, then we even have equality. Now assume E[Z] > 0. By concavity of h, there exists a 
b e R such that for all z > 0, 

<h(E[Z]) + (z-E[Z])b. 
Since h is nondecreasing, we have b > and thus 

E|Y/i(Z)] < E[y (h(E[Z]) + (Z -E[Z])b)] < E[Y]h{E[Z]). □ 

Let l\,h € 5 and let X K be the truncated process with initial state Xq = [ts^x, tn 2 \). 
Writing (f2T2T|) in the form 

+ £ / a t - s (k,l)dM?(l), 
i&s Ja 

and recalling that the M^(l) are orthogonal martingales, we get that the random variables X^ s (k) and X^ s (l) 
are uncorrelated for k,l £ S. Hence aX^ s (l) and s (/) are uncorrelated. Note that x 1— > hx{x) (defined in 
(|2.20[) ) is concave and monotone increasing for x < K. Hence we get by Lemma \2. 61 that 

E[(Mf(0> t ] < t hK(E[X^(l)])E[aXl jtS (l)]ds. (2.25) 
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Denote by 

Mf =^Mf(fe) = (Xf,l>, i = l,2, (2.26) 

the total mass process of .X^ . Then (|2.25p implies 

Var[M^]=E[(Mf>J 

= ^E[(Mf(0) ( ] 

ies (2.27) 



2 

- S / X! ^((a«J0,o)(0) ( a ' a^3-j,o)(0 ds. 

.•-WO 



j=i" u les 

Using again the concavity of Iik and Jensen's inequality (for the probability measure 1 1— ► (a • a s )(l, k)) we get 
Var[M£] < £ / , °°^X 3 _,-o(fc)^(^(a-a s )a,fc) (a a X;,o)(o) ds. (2.28) 
Now recall that the initial states are X^o = and that a s = (a s (a + a T )a^). Hence we obtain 

2 poo 

Var[M^]<V/ h K (a s (h,h))ds <8Klog(K)G loe (hJ 2 ). (2.29) 

3=1 - 70 

Hence M/^ and M% are (orthogonal) L 2 -bounded martingales and thus converge almost surely and in I? to 
some random variables and Mf^ with E[Mf] = 1 and VarfM^J < 8K \og(K) G\ og (h,l 2 )- 

2.4 Proof of Theorem [j] 

Clearly, the product ■ is a uniformly integrable martingale, hence 

E [M^M&J = E [M*M*] = 1. (2.30) 

Write 

Mi := M i;00 1{ Tk<oo} . 

bincc on t k = oo, in order to show Theorem [IJ it is enough to show that EfMff^Mj^ l {r/f=0 o}] > 0. 

To this end we compute 

EfM^M*^ l {rjr=oo} ] = E[M^M£J - EfM^M^ 2 ] - E[Mf MfJ + E[Mf M 

By (|2.30|) . it is thus enough to show 



EtM^Mf ] < 1 and E[M* M£J < 1. (2.31) 

Although Miqq and ^ are uncorrelated, this is not true for M-^ and (at least we cannot show this). 

Hence we have to use a slightly more subtle argument. We employ the Cauchy-Schwarz inequality to estimate 

E[M^Mf ] 2 < E^) 2 ] E[(Mf ) 2 ] 

= (1 + Var[M*y)E[(Mf ) 2 ] (2 32) 

3K x 2 
IT 



< i + 8in og (i<:)Giog(M2; 



Recall that X > 2 is fixed. Now choose l\ and ?2 such that Giog^i,^) becomes so small that the right hand 
side in |2~32| is bounded by \. Similarly, we get E[Mf M^f < \. 

This shows (j2.31[) and thus completes the proof of Theorem [TJ □ 
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3 Non-coexistence of types, Proof of Theorem [2] 



The strategy of proof is described in the following two steps. 

Step 1. Replace the process X by the approximate process X K constructed in Section [2T2l If X would have 
coexistence of types, then so would X K (for some large K). 

Step 2. Since M^ t = (X^ t , 1), t > 0, is a convergent martingale with bounded jump size, also the conditional 
quadratic variation process (M^. ) converges. We derive a lower bound for (M/^.) and show that due to the 
recurrence of A, this lower bound would diverge with positive probability if X K had coexistence of types. Hence 
there can be no coexistence of types for X K and thus neither for X. 

3.1 Step 1: The approximate process 

Assume that for X, coexistence of types is possible. We want to lead this assumption to a contradiction. 

Recall that Mn := (A;. t , 1), t > 0, i = 1, 2, are the total mass processes. Coexistence of types means that there 
exists a deterministic initial state Xq such that Mi t t < oo, i = 1, 2, and such that 

lim Mi t M 2 t > with positive probability. 

t — >oo ' ' 

(Recall from the discussion prior to Theorem [T] that the total mass processes are nonnegative martingales and 
are hence convergent.) That is, there exists a 5 > and a To such that 

P[M M M 2)t > S for all t > T ] > 65. 

Imposing an optional sampling argument, we see that (Mi t t)t£[Q,T ] is bounded away from with probability 1 
if Mi t T > S. Hence, with a possibly different 5 > 0, we get 

P[M M M 2 , t > 5 for all t > 0] > 68. (3.1) 

Now define 

K :=l(M lfi + M 2fi ). 

Let X K denote the truncated process defined in Section l2~2l with X$ = Xq. Recall that (M^ t )t>o, is the total 
mass process of Xf- and that it is a martingale. Hence we have 

P[A]>1-S, (3.2) 

where 

A := {A/5 + Mg t < K/2 for all t > 0}. (3.3) 

Also 

P [M M + M 2 ,t < K/2 for all * > 0] > 1 - S. 
We can couple X K and X such that both processes coincide on A. Hence for 

B := {Mg t Mf; t > S for all t > 0}, 

we have 

P[B] > P[4nB] > P[M u M 2 , t >S for all t > 0] - P[A C ] > 55. (3.4) 
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3.2 Step 2: The lower bound for the conditional quadratic variation 

Denote by {{M\.) t )t>o the conditional quadratic variation process of (M^ t ) t >Q. Since is a martingale whose 
jumps are bounded (by K), convergence of implies almost sure convergence of — > < oo as 

t — > oo. Hence for any e > and S > 0, there exists a To such that 

P [( M V> t -( M V>T > £ ]< S fora11 *> T o- ( 3 - 5 ) 



The aim of this section is to show that (|3.4p leads to a contradiction to (|3.5p which shows that the assumption 
that for X coexistence of types would be possible was wrong. 

Fix S > 0. We choose an appropriate e > for a contradiction via the following procedure. Recall (|1.19p and 
define 

a*(k) := 1- a(k,k) and a* := inf a*(k) > 0. (3.6) 



Recall c from (|1.20p and define 



(3-7) 



Let 

2 a: 2 

^=2+—. (3.8) 

Furthermore, define 

£: =tTT (3 ' 9) 

Note that 

K 2 

e - c'S 2 = -Re and — < 6. (3.10) 

Recall from (12.21)) that M^(k), k £ S, are orthogonal martingales and that = YlkeS-^i s )- Hence, by 
(l2~24l) . we have for t > T , 

{M*) t - {M«) To = E « M i*( fc ))t - < M i K ( fc ))T ) > o E / + oX* (fc)X^(fc)) 

fees kes jToATK 

where is defined in (|2.23[) . Define 

1 rt 

Z t = - V / Xf s (A;)^Xf s (fc) + X%.(k)AX?„(k) ds for * > T . (3.11) 

2 fees Jt ° 

Note that t k = oo on A. Hence for all t > T , 

(Mg) t - (Af£) To > Z t on A. (3.12) 

Lemma 3.1 Tor any k±,k 2 € S, k\ ^ k 2 , let N t (ki,k 2 ) be given by 

Nt(k u k 2 ) = E f a ts(ki, l')at-s(k2, 1) (X« s (l')dM* s (l) + X« s (l)dM*(l')) for t > T . 

Ml' 

Then for t > To, on A we have 

XZMXZtfa) = at -T ^ To (fc 1 )a t _ To ^ To (fc 2 ) 

-E / at-s(kul)a t . s (k 2 ,l) (X* s (l)AXK s (l) + XK s (l)AX* s (l)) ds + N t (ki , k 2 ). 
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Proof. This is an immediate consequence of (I2.21[) . the fact that 

Xg t (k)X$; t (k) = for all ke S, if t < t k , 
and since t k — oo on A. □ 

Lemma 3.2 On the event A, the following decomposition holds: 

Z t = Z\ 1] - Z\ 2) + Zfl + Zfl for t > T , (3.13) 



where 



Z ^ = E / a -*>*Sr G (0 (Aas-T X« To )(l) +a*(l)a s - To X* To (l)) ds 
Z t 2) = /ZY,f I a s - r (l,k)[(Aa s - r (;k))(l)+a*(l)a s -r(l,k)] 



ies kes jT ° jTo (3-14) 
(x5(fc)aXf r (fc) + aXf r (fc)X£.(fe)) ds, 



Z T1 = E S C> ^ E / / a — C' fc ) (°-r*£i,r(0 - "-rtf'. *)*3* i,r(*)) dM?(k) ds. 

Proof. By definition, 

Z * = ^EE f X?{k)AX*_ ita (k)d s (3.15) 

= \ EE f X^{k)Y,Hk,l)-t {k=l} )X^_ itS {l)ds 
i=ikes T ° les 

= t a(ki,k 2 )X« s (k 1 )X«Jk 2 )ds, 

k 1 ,k 2 es Jfo 

where the last inequality follows since X^ < t (k)X 2 K t (k) — for all k S 5, t > on A. Now by Lemma [3.11 we 
have 



Z t = E / a ( k i, k 2)a s -T X^ Tg (k 1 )a s ^ To X 2 K To (k 2 )ds 

k lt k 2 £S JTq 
k\^k 2 

- E / a(fci,/c 2 )E / a s-r{kl,l)a s - r (k 2 ,l) 



(3.16) 



x (X5(0^X 2 K r (0+X 2 K r (/)^Xf r (/)) drds 
^ / a(fci, k 2 )N s (ki, k 2 ) ds. 

JTn 



hi ,k 2 £S 

k\=£k 2 



Let us consider the first term on the right hand side of the above equation. We easily get that 

Y] / a(fcl,&2)a s _To^T (M a s-To^T ( fc 2)ds 



k 1 ,k 2 £S 
k\^k 2 



(3.17) 



fci£S 
7 (1) 
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and we are done with the first term in the decomposition. Similarly, for the second term we have 
^ / a(fci,fc 2 )y^ / a s - r (ki,l)a s _ r (k 2 ,l) 

x (x£ r (l)AXg r (l) + Xg r (l)AX£ r (l)) drds 



k 1 ,k 2 es J T leS jT ° 

fc l7 £fc 2 



= / O.-r(*l,0[-Al»-r(-,0(*l)+ o *(*l) a «-r(Ali0] 

x (xf r (o>ur^ r (o + ^ r (OAyifr(0) drds 



'l 2) 

For the third term, we have: 



J2 f a(k 1 ,k 2 )N s {k 1 ,k 2 )ds (3.18) 

1 5^ fe 2 

= J2 [ Hki,k 2 )J2 f ^(ln.Oa,-,^,!) (1^(0^(0 + ^(1)^(1')) els 

fc i?^ fc 2 i^i' 

^ a(k u k 2 ) f f |^a s _ r (fc 2 ,Z)(a s _ r X 1 ^(fe 1 )-a s _ r (A; 1 ,Z)X^(0)dM 2 i f r .(0 

fci,fc 2 es •'To ''To [ ;eS 

fel ^&2 

J] a._r(fcl,0(a.-r*v(*2) " O.-r(*2,0-yv(n)^v(0 | ds 

J'es J 
53 B(fci,te) / / |53a s _ r (fc 2 ,Z)(a s _ r X 1 ^(fe 1 )- as _ r (A; 1 ,Z)X^(0)dM 2 ^(0 

i>fcaes ■'To ■'To [ ; GiS 



53a s _ r (fc 2 ,Z)(a s _ r X|f r (A ;i )-a s _ r (fe l! Z)X 2 x r (0)dM 1 i f r (Z) L 



_ 7 (3) , 7 (3) 

where the equality before the last follows by the symmetry of a. □ 
Lemma 3.3 Recall G* from (|1.15|) . For aZZ i > Tq, on A we have 

z (2) < G *2{t~T ) ^ 

Proof. Recall (|1.13p and (|1.14p and note that (with / the unit matrix) 

* - 1 

a T r Aa r dr = -{a 2t - I). (3.19) 
o ^ 
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Hence (using that a* (I) < 1) 



Z[ 2) < V f f [a T s _ T Aa s - r + a 2{s _ r) ] (k, k) ds (X^ r (k)aX^ r (k) + aX? r (k)X* r {k)) dr 

/"* / 1 /"* \ (3 20) 

= ^2 J \^2{o,2(t-r){k,k) - l) + J a 2 ( s - r )(k,k)dsj 



fees' 

[Xg r (k)aX% r (k) + aX^ r {k)X^ r {k)] dr. 

Then use the fact that a t (k, k) - 1 < and aX^ r (k)X^_ . r (k) = AX^ r (k)X^_ tr {k) on A to get that on A, zf ] 
is bounded above by 



\ E / (AX* T (k)X* r {k) + X* r (k)AX* r (k)) dr 



/T ° fees 



Next we will handle Z; . 



Lemma 3.4 i<br all t > and /or i = 1,2, we have 



kes JT ° 

where for t sufficiently large 

tK 



Z ti = E / fci(M,r)dM£(A), (3.21) 



|^(fc,t,r)| < G* (i _ To) M 3 ^ iir /or a/Z r e [T ,t]. 

Proof. First, by the stochastic Fubini Theorem (see, e.g., Theorem 64 in [H]), we can change the order of 
integration in order to get 

hi(k,t,r)= S (M') / a a - r (l,k)(a s ^ r XK_ ir (l') - a s - r (l',k)X£L i>r (k)) ds 



t 

n 

s — r 
t 



les 
=: h+h- 
For the first integral, using (|3.19|) we get 



a; _^Aa s - r X s _ i r (k) — a s _ r Aa s - r (k, k)X 3 _ i r (k) ds 
J Yj a*(l)a s - r {l, k)[a s - r XK_ i<r (l) - a.- r (l, k)X^_ %r {k)] ds 



< ^a 2it _ r) Xl t Jk)-^a 2(t _ r) (k,k)X^ h M = h < \M$_ i>r . 



For the second integral, we obtain similarly 



< h < f E as ~rV, k)a s - r XgL iir (l) ds 

Jr t^a 



zeS 

t 

K 



a s - r X 3 _i r {k) ds 



1~ 1 
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Combining the estimates for I\ and I2, we get 

Mk,t,r)\ < \{Gt {t _ To) +l)M^ r < G* 2{t _ To) M 3 K _ hr , 

and we get the bound for hi(k,t,r) for t sufficiently large. □ 
We have to introduce some notation and define a number of additional constants. Define the event 



C = \ E X *T (k) > \m* To for i = 1,2 I . 

[\k\<L J 



(3.22) 



Since Mf- To < 00, i = 1, 2, there exists an L > such that 

P[C]>l-<5. (3.23) 

Lemma 3.5 There exists a T\ > To suc/i i/iai /or a/Z i > T\, 

Z[ 1] > c> Gl {t _ To) M« T M* To on CnA. 

Proof. To simplify notation let z\ 1] = Z t (M) + Z t (1 ' 2) , where 

: = E / a.-ib^7b(iM(a.-ib^T 11 )(0'fc 1 

zes •' T ° 

We start with showing that Z^' 1 ^ > 0. To this end we compute using (|3.19p 

Z t (M) = f X«T Q) [a T s _ To Aa s -T ]X« T() {l)d S 
les T ° 

= \ E^ o (0(S 2 (t-To) - ^)^T (0 

zes 

= ^E^o(0«2(t-T )^T (0 > On A, 

zeS 

since (fc)Xf To (fc) = on A for all k E S. 

The bound for Z^ 1 ' 2 ^ follows similarly to Dawson-Perkins p. 1109]. Recall c from (|1 .20|) . For k,l € S such 
that \k\, \l\ < L, let T(k,l) be large enough such that 

^M>| for all t>T(k,l). 

Define 

Ti = T + max T(fc. f) . (3.24) 

|fc|,|Z|<L 

Then for any t > T\, we have that on C 

> y G 2(t-T )(k,l)X^ To (k)X^ To (l) 

k,ies 

^ yf E ] min G 2(t _ To) (fc,Z) 



a c . 



- ~Jq M 1,T M 2.T a G *2 ( t - T ) ) 
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where the last inequality follows by (|3.22p and (|3.24|) . Recalling (|3.7|) we get that on C 

zj 1 ' 2 > c' G2( t _ To -)M^ To M^ To . 
Since z£ > on A, this finishes the proof of Lemma 



From Lemma 13731 and Lemma 13751 we get that oniflC, Z t is bounded below by 

G 



2(t-To) Zt + jj( 3 ) 



Z t > c G* 2 ^ t _ To ^M l ^ To M 2 To 
where zf ] = zf} + zfj. Let a := 1/(1 + G* 2(t _ Tg) /2). From ([3~2~5]) we get that on A n C 

Then (recall from (|3~T2"j) that Z t < (M£) t - (M£) on A) 

P[(M 1 ^> t -(M 1 J f.) To < £ ,AnC] 

= P [Z t < e, (M*) t - (m£) <MnC] 



< P 



a Z t (3) < £ - ac' G* 2{t _ To) MK To M* To , (M 1 K ) t - (m£> < e, A n C 



We assume that i is sufficiently large so that G 



2{t-T ) 



> 2, hence 



1 < aG 



2G 



2(t-T ) 



2(t-T D ) 2 + G. 



< 2. 



2(t-T ) 



By ([3728]) . pTi]) and (f3TT0|> . we get 



aZf ] < e -ac'G* 2{t _ To) M^ To M 2 K To , (M*) t - <Af£) < e, A n C 



< P 

< P 
= P 



aZ\ 3) <e-ac'G* 2(t _ To) 5\ (Af£> t - (M^ < e, AnC 



+ 1-5(5 



Z t (3) < -ife, (Af£) t - (M£) T < e, A n C 



< i?" 2 e~ 2 a 2 E 



( 4 3) ) 2 1 { < M * > t - < ) ^ <e} *A 



1 - 5(5 
1 - 5<5 

1 - 5(5. 



Using Lemma 13.41 this inequality can be continued by 



<i?- 2 e- 1 (if/2) 2 (aG; (t _ To) ) 2 + l-5J 



K 2 
< 1 - 4(5. 

Combining (j3~27j) . 1(3725]) . (j3~2|) and (|3~23| . we get 

P[<M*> t - (M£) n <e]<P [(M 1 K ) t - (M*) To <e,Anc]+ P[A C ] + P[C C ] 

< 1-4.8 + 6 + 5 = 1-25. 
This is a contradiction to (|3.5|) and hence finishes the proof of Theorem [2j 
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